Coupled thermo-mechanics of single-wall carbon nanotubes 
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Abstract 

The temperature-dependent transverse mechanical properties of single- walled nanotubes are stud- 
ied using a molecular mechanics approach. The stretching and bond angle force constants describ- 
ing the mechanical behaviour of the sp 1 bonds are resolved in the temperature range between K 
and 1600 K, allowing to identify a temperature dependence of the nanotubes wall thickness. We 
observe a decrease of the stiffness properties (axial and shear Young's modulus) with increasing 
temperatures, and an augmentation of the transverse Poisson's ratio, with magnitudes depending 
on the chirality of the nanotube. Our closed-form predictions compare well with existing Molecular 
Dynamics simulations. 
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Carbon nanotubes-based composites, both based on metal/ceramics or polymer ma- 
trix undergo substantial thermal loading during both the manufacturing phase and their 
operational life. This specific aspect of their multifunctional behaviour has generated a sub- 
stantial amount of activity in the identification of the thermal properties of nanotubes, with 
particular emphasis on the thermal conductivity (3-7]. The dependence of the Young's mod- 
ulus and tensile strength of single- walled carbon nanotubes (SWCNTs) has been investigated 



using various MD techniques 8Mll| and hybrid atomistic-Finite Element techniques 12l. 1 13 [ - 
The importance of the interplay between the nanotube and the polymer matrix during differ- 
ent environmental temperatures has been highlighted by Wei from Monte Carlo/Molecular 
Dynamics (MC/MD) simulations on a periodic assembly of a PE-(19,0) CNT nanocompos- 
ite jl^). The use of simulations to predict the temperature dependence of the nanotubes 
mechanical properties is a necessity, because of the difficulties involved in obtaining direct 
experimental measurements in different thermal environments. However, all the predictive 
tools so far available to simulate and design the stiffness of CNT-based composites are not 
based on analytical formulas that would facilitate the work of the material scientist and 
engineer. The existing atomistic-continuum models describing the elastic properties of car- 



bon nanotubes in closed form 
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16| have their bonds force constant formulated only for 
room temperature. In this work we propose a compact formulation that express the re- 
lation between the sp 2 C-C bond force constants and the surrounding temperature of the 
nanotube, providing a functional description of the coupled thermal and mechanical linear 
elastic properties of SWCNTs. 

In the following description of the model, we will use the suffix T to indicate the 
temperature-dependent force constants, while the suffix ° will indicate a force constant 
term referred to the room temperature (T = 300 K ) . Adopting in part the Universal Force 
Field (UFF) nomenclature 17j, the stretching force constant related to the calculation of the 
harmonic harmonic potential energy associated to the stretching deformation is expressed 
as: 



664.12 



Z*Z 
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Where Z\ and Z*j are the effective electronic charges (1.914 electron units). At room 
temperature, the equilibrium length r°j between the atoms / and J is defined as 0.142 nm, 
with the stretching constant ([T]) assuming the value of 6.52 10~ 7 N nm' 1 , which is consistent 



with the force constant used in the AMBER model (is| . The bending force constant kJ JK 
related to the harmonic potential energy associated to the bond bending angle is the second 
partial derivative of the bending energy functional 



k ui< = (~}JrJ = ~ Z J^ K riJ rjK { 3r v r JK [ l ~ (cos9) 2 ] - r 2 IK cos#} (2) 

Where the distance between the atoms / and K is expressed as tjk = t\j + r 2 K — 
ZrjjrjK cos 6. At room temperature, the bond angle is equal to 9 = 120 C°. The term 
(3 is equal to 332.6/rjj/rjk to obtain at room temperature a value of kf JK consistent with 
the one provided by the AMBER force, and used by several authors in atomistic-continuum 
mechanics models (8.76 10~ 10 N nm/rad 2 ). The torsional potential energy between bonds 
I J and KL is approximated as the summation of three cosine Fourier terms based on oj (the 
angle between the IL axis and the UK plane): 
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At room temperature the constant Kj JKL is assumed equal as 2.78 \0~ lo N/ nm/rad 2 , 
consistent with the value used by several Authors fiol . [^j. Due to the very low sensitivity 
of the torsional constant versus the temperature, in our model we consider the term Kf JKL 
as constant with the different thermal environments, and will be indicated as Kj JKL in the 
rest of the work. Following Chen et al jl3|, the bond angle 9 is also assumed constant at 
the room temperature value (9q = 2tt/3). The thermal dependence of the force constants 
(JTJ and (j2J) is provided by the thermal variation of the bond length through an equivalent 
coefficient of thermal expansion (CTE) a: 

rjj^rjjil + aAT) (4) 

Where AT = T — Tq. In this work, we consider as CTE for the bond the coefficient 
of thermal expansion of suspended graphene sheets calculated using a non-uniform Green's 
function approach by Jiang et al |2l|. At room temperature, the coefficient of thermal 
expansion is equal to -6 X lO -6 /^" 1 , 14 % lower than the experimental value measured 



by Bao et al 22|. The CTE assumes a value of -3.2 X lO -6 -^ -1 at 20 K, to decrease 
to a minimum of -1.25 X lO -6 -^" 1 at ~ 80 K. After this temperature, the coefficient of 
thermal expansion increases monotonically, although it remains negative up to 650 K. At 
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T = 1550 K, a brings a value of 5.0 X lO -6 ^ -1 . Assuming rjj = rjx, using (j3J) in ([I]) and 
(|2J), it is possible to express the harmonic potentials related to the stretching and bending 
energy respectively with the variation of the environmental temperature. 

The equivalent mechanical behaviour of the C-C bond is represented at this stage by 
equating the harmonic potentials with the correspondent strain energies deformations under 
axial, bending and torsion of a structural beam element 19]: 
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In (J5J), E and G are the equivalent Young's and shear modulus of the material representing 
the C-C bond, I = 7r<i 4 /32 is the polar inertia moment of the bond beam (considered 
having a circular cross-section of diameter d, equal to the thickness). Differently from other 
approaches used in open literature 23|, we adopt in ([5]) a Timoshenko beam having deep 
shear cross-deformation behaviour as representative beam model, to consider more realistic 
distributions of thickness and equilibrium length existing in graphene and nanotubes 



The effect due to the shear-induced cross section deformation is provided by the constant 
$ = 12EI/G/A s /rjj 
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where A s = nd /A/F s is the reduced cross section of the beam 



by the shear correction term F s depending on the Poisson's ratio v of the equivalent C-C 



bond material 



29J: 
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The equivalence between the harmonic potentials and the beam strain energies in 
leads to the following set of equations: 
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Inserting ([6]) into the definition of the shear constant $, and solving for E in (J2J), we 
obtain the following nonlinear expression in d and v for a specific temperature T: 
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We impose the additional condition that the equivalent material of the C-C bond beam 
behaves as isotropic, i.e., with no directional preference in its mechanical response. The 
isotropic condition G = E/2/ (1 + u), together with relation (jHJ) constitutes a system of 
nonlinear equations which can be solved with traditional methods, such as the Marquardt 
algorithm. The solution of the system yields a unique value of the thickness d for a given 
temperature T. At room temperature, the thickness identified using this approach yields 



a value of 0.084 nm, equal to the one found for single wall and nanotube bundles in |25|. 
We observe a low sensitivity of the thickness corresponding to the isotropic condition versus 
the temperature, with a minimum value of 0.0835 nm at 1600 K, and a maximum value of 
0.0841 nm for T = 470 K. The temperature dependence of the thickness can be described 
using a polynomial curve of the 9 th order, as shown in Table HI It is interesting to notice that 
the C-C bond behaves mechanically as if made of a quasi- zero Poisson's ratio behaviour, like 
in natural-occurring cork ^30 1 . The 0.084 nm value at room temperature is also consistent 



with the 0.084 nm found by Kudin et al 



3lfl , and 0.074 nm identified by Tu and Ou-Yang 



32| . From simulations using the MM3 potential, Batra and Sears have found an equivalent 
thickness of 0.043 nm for different chiral configurations 33|, although a value o: 



0.1 nm 



was observed in nanotubes undergoing bending and breathing vibration modes [34j. Zhang 
and Shen identify thickness of 0.088 nm and 0.087 nm for (17,0) and (21,0) nanotubes 
respectively, while for (10,10) and (12,12) the thickness calculated is 0.067 nm 10j. The 



thickness values calculated with our method do compare well with the 0.080 



and Cao 



although they are higher than the 0.066 nm of Yakobson et al 



nm of Vodenitcharova and Zhang 



37|. 



nm of Chen 



36|, and 0.60 



Classical mechanistic theories for the elastic properties of single-wall carbon nanotubes 
consider the strain energy associated to the deformation of sp 2 bonds in terms of stretching 
and bond angle constants [C p and Cq respectively 15j,LL6|). For the bond angle constant (mo- 
ment), we consider that for linear elasticity the lattice is dominated by hinging deformation 
38|. In cellular structures under general flexural/axial deformations, the hinging constant 
can be calculated as Cq = El/q, where q is the portion of the C-C bond length undergoing 



hinging, while the segment rjj — q deforms as a rigid body 



39|, see also Figured]). Using 



the relations ([7]), the bending angle constant can be rewritten as Cg = kjjd 2 rjj/16/q. Con- 
sidering a ratio rjj/q = 5, inserting the numerical values of the thickness and equilibrium 
length at room temperature we obtain a value of 1.44 nN nm rad~ 2 , well in line with the 1.42 



nN nmrad 2 used in 
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16| . The stretching constant is calculated as C p = k p kjj, where 



hp = 1.128 to converge to the 0.36 TPanm 1 of the in-plane graphitic surface modulus (chi- 



ral index n — > oo] 
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40|. At room temperature, the stretching constant identified in this 
model is equal to 735 nN nm -1 , which compares well with the 742 nN nm -1 used by Shen 
and Li 16|. Cadelano et al f4fy have identified a value of the surface modulus of graphene- 
type systems equal to 0.312 TPanm~ l at zero temperature damped dynamics using tight 



binding calculations and continuum elasticity. Using our approach, we observe a value of 
0.332 TPanm -1 at zero temperature, 6 % higher. For both the bending and stretching 
constants, we observe a very low dependence versus the environmental temperature between 
370 K < T < 570 K. The stretching constant then decreases with increasing temperatures 
with an approximate rate of 0.015 nN nm~ x K~ x , a behaviour similar to the one of the Cg 
constant, having the latter a decrease rate of 4.9 10~ 5 nN nmrad~ 2 K~ l between 600 K and 
1600 K. At low temperatures, both the force and bending constants decrease from values 
equal to 99.2 % of the maximum Cg and C p at 1 K to a minimum at 70 K corresponding 
to the 98.6 % of the maximum of the constants at 450 K. Similarly to the C-C bond thick- 
ness, the force stretching and bending constants can be expressed in polynomial terms, as 
illustrated in Table [H 

The temperature-dependent constants C p and Cg can now be used in atomistic-continuum 
approaches illustrated in open literature to model in a closed-form solution the mechanical 
properties of SWCNTs under different temperature conditions. In this work we have used 



the formulation proposed by Shen and Li 



16| to simulate the axial and shear stiffness, as 



well as the Poisson's ratio of the nanotubes with armchair and zigzag geometry. We compare 



the surface Young's modulus Y{ x calculated in analytical form against 
through MD simulations using a REBO potential by Zhang and Shen 



the results obtained 
101 (Figure ED. Our 



predictions arc well in line with the 0.352 TPanm 1 of Tu and Ou-Yang |32| and Pantano et 



al 42], although 3.5 % higher than Li and Chou 19] and Wang and co- Authors [43]. In terms 
of longitudinal modulus Yj ± (defined as Y^/d), our predictions tend to overestimate the MD 



simulations in 



10j. For example, our axial moduli for (17,0) nanotubes are 9 % higher than 



the molecular dynamics simulations at T = 300 K, while a at a temperature of 1000 K our 
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overestimate is around 11 %. The effective Young's modulus Yn = Yj / (i?/2)(where R is 



the diameter of the nanotube) has been also derived by Hsieh et al |ll| from the vibration 
amplitude of (n, 0) clamped-free nanotubes simulated with a Tersoff-Brenner potential at 
temperatures ranging between K and 2000 K. The MD simulations in {h| are least-squares 
fitted to an asymptotic value of the graphene Young's modulus equal to 1.2 TPa found by 



44 



45] at a not specified temperature, 20 % higher than the experimental value found 



w Lee et al for a thickness of 0.34 nm 



46] . A general comparison of the results from 



ll| against our simulations is presented in Figure [21 showing a general agreement between 



trends related to the nanotube radius at T = 1100 K, and within the temperature range 50 
K - 1100 K. We observe a good convergence between MD and our simulations especially for 
CNT radius lower than 0.5 nm, and excellent agreement (percentage error less than 4 %) for 
(7,0) SWCNTs over the temperature range considered. Hsieh and co-Authors record a steep 
decrease in terms of axial stiffness above 1100 K, corresponding to an abrupt change of the 
standard deviation associated to the dynamic tip displacements of the tube We observe 
also an overestimate for the shear modulus G\ 2 flO|. with discrepancies around 29 % and 
28 % at 1000 K and 1200 K respectively. However, our predictions show a decrease of the 



longitudinal shear modulus versus the temperature, opposite from the results derived in [lOj. 
The transverse Poisson's ratio v 2 \ increases with the temperature (Figure HJ), although the 
augmentation is limited on average to a 3 % increase between T = 300 K and T = 1200 K 
for the different chirality and tube diameters. The low sensitivity of the Poisson's ratio over 
the temperature is due to the weak dependency versus the temperature itself of the term 
Cprfj/Ce, which is contained both in the numerator and denominator of the Poisson's ratio 



expression 



16|. Chen et al. 13J observe also an increase of u 2 i with temperature (31 % of 



increase for zizgag tubes of 1.2 nm diameter between room temperature and T = 1200 K). 
The discrepancy between theirs and our predictions is due mainly to the coefficient of thermal 
expansion used for the C-C bond, being positive in their case for T > 298 K, while our model 
has a positive CTE for temperatures higher than 650 K |2l|. We observe that Zhang and 



Shen predict a transverse Poisson's ratio for a (17,0) nanotube at T = 300 K equal to 0.17, 



in good agreement with our value of 0.16 



lOj. We note that for graphene- type systems, 



our model predicts at zero temperature a Poisson's ratio of 0.12, significantly lower than 



the 0.31 identified in 



41]. However, our range of Poisson's ratios at room temperature for 



zigzag tubes (Figure 2]) is in line with the 0.12 - 0.19 values predicted by Sanchez- Portal et 
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al 



47j, the 0. 



Zhou et al 



49(]. 



5 of Kudin and co- Authors |31fl, 0.21 of Sears and Batra [48], and 0.24 from 



In summary, our model allows to predict in a compact form the overall elastic trans- 
verse properties of single wall nanotubes over a large temperature range. The bending and 
stretching constants identified with our approach capture the intrinsic hinging/stretching 
behaviour of sp 2 C-C bonds in nanotubes, and relate the stiffness and Poisson's ratio to the 
ambient temperature of the nanotube, with 2 % average decrease of the Young's modulus 
between the ambient temperature and 1200 K for a given tube chirality. 
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Quantity 


Pi 


P2 


P3 


Pa 


P5 


P6 


Pi 


P8 


P9 


Pio 


d X 10" 5 [nm] 


-2.256 


2.571 


12.7 


-13.6 


-17.78 


13.14 


14.81 


14.13 


-23.66 


8401 


C e X 10 12 [N nm rad~ 2 ] 


-1.84 


2.181 


9.427 


-9.675 


-14.92 


12.9 


12.47 


-11.97 


-19.8 


1435 


C p X 10 10 [Nnm- 1 ] 


-5.731 


6.702 


2873 


-29.83 


-45.55 


35.13 


37.83 


-36.98 


-59.98 


7340 



Table I: Polynomial interpolations for the thickness, stretching and hinging force constants for the 
sp 2 C-C bond lattice versus the temperature. The polynomials are of the type y = Xm=i PnT w ~ n , 
where f = (t - t) /a (T), with T = 800 A" and standard deviation a (T) = 466 K. All fittings 
have an average R 2 of 0.991, with 95 % of confidence interval. 




Figure 1: The variation of the stretching constant C p and hinging constant Ce over the temperature. 
The hinging constant is calculated considering the change in angular deformation 9 — 6' = AO = 
Mj dq = Mq/EI. Under hinging, only the initial portion q — rj j/5 of the bond length deforms 
by rotation, while the rest of the bond deforms as a rigid body [39fl. The hinging constant is then 
defined as C e = M/A6 = El/q. 
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Figure 2: Comparison between the surface Young's modulus of different zigzag and armchair SWC- 



NTs and the MD simulations from 
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Figure 3: The effective Young's modulus modulus for (n, 0) SWCNTs. The results from the proposed 



method (continuous line) are compared with the ones from Hsieh and co- Authors 



11| (circular dots). 
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Figure 4: Map of the transverse Poisson's ratio z/21 for (n, 0) tubes versus the the chiral index n 
and temperature. 
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